Mon. Not. R. Astron. Soc. 000, (2000) Printed 1 February 2008 (MN I*TeX style file vl.4) 



Using galaxy redshift surveys to detect 
gravitationally-lensed quasars 

Daniel J. Mortlock^'^'^^ and Rachel L. Webster^^ 

^School of Physics, The University of Melbourne, Parkville, Victoria 3052, Australia 
'^Astrophysics Group, Cavendish Laboratory, Madingley Road, Cambridge CB3 OHE, U.K. 
^Institute of Astronomy, Madingley Road, Cambridge CBS OHA, U.K. 



o 
o 
o 



Accepted. Received; in original form 2000 April 19 



ABSTRACT 

Gravitationally-lensed quasars can be discovered as a by-product of galaxy redshift 
surveys. Lenses discovered spectroscopically in this way should require less observa- 
tional effort per event than those found in dedicated lens surveys. Further, the lens 
galaxies should be relatively nearby, facilitating a number of detailed observations 
that are impossible for the more common high-redshift lenses. This is epitomised by 
the wide range of results that have been obtained from Q 2237-1-0305, which was dis- 
covered as part of the Center for Astrophysics redshift survey, and remains the only 
quasar lens discovered in this way. The likelihood of this survey yielding a lens is cal- 
culated to be ~ 0.03, which is an order of magnitude larger than previous estimates 
due to two effects. Firstly, the quasar images themselves increase the observed flux 
of the lens, so that lens galaxies up to a magnitude fainter than the nominal survey 
limit must be included in the calculation. Secondly, it is possible for lensed quasars 
with extremely faint deflectors to enter the survey due to the extended morphology of 
the multiple images. Extrapolating these results to future surveys, the 2 degree Field 
galaxy redshift survey should contain between 10 and 50 lenses and the Sloan Digital 
Sky Survey should yield between 50 and 300 lenses, depending on the cosmological 
model and the observing conditions. 
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_C3 1 INTRODUCTION particularly useful, but for rather different reasons. The best 

„ , , , , , , . illustration of this is Q 2237-1-0305 (Huchra et al. 1985), a 

Gravitationally-lensed quasars are extremely valuable tor i i . i , , i . r. r i . i i 

, r 1 quadruply-imaged quasar at a redsniit oi 1.69, seen through 

a number ot reasons. The frequency oi lensme events is a : , > i r • , , , j i -^t r r, n ^ r^n 

, , the bulge oi a spiral galaxy at a redshift of 0.04. ihe prox- 

sensitive probe of the cosmological constant (e.g. Turner . r j_i i i • i j_ -j, • i j_ j_ i^i • tt 

, \ , , , r , imity ot the lens galaxy is such that it is detectable m 11 1 

1990: Kochanek 1995, 1996a) and the density ot the uiii- , , T„r>r>\ i i- n j i j- 

' (Barnes et al. 1999), has an optically measured velocity dis- 
verse (Kochanek 1995; Mortlock & Webster 2000a). Indi- . , , irvr,r,\ i i i i i 

\ , persion (loltz et al. 1992), and can be resolved on scales 
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„ , , r, / ^, , , n TT . -.^^^ ^7 01 ^ 100 pc by the nubble Space lelescope. ihese results 
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Turnshek & Khersonsky 1998; Keeton, Kochanek & Falco 
1998), the mass of galactic halo objects (e.g. Schmidt & 



lensing constraints to provide a wealth of information about 

. , , „ ^ , , both the galaxy and the lensed quasar. The macroscopic 

Wambsganss 1998; Wyithe, Webster & Turner 2000b), ex- ^- c ^1 ■ „ ^ j ^ ■ ^- 

" / o properties ot the images allow accurate determinations of 

tinction in the lens galaxy (e.g. Malhotra, Rhoads & Turner i- ^ -t, j_- 1 ^ i- 1 j_ j_- r j_i 1 

„ , , , \ „ , , , , , „ „ , , the mass distribution and mass-to-light ratio ot the galaxy 

1997; Falco et al. 1999), and Hubble's constant (e.g. Refsdal ^ 5 ,y 



1964; Grogin & Narayan 1996; Kundic et al. 1997). Further, 



(e.g. Yee 1988; Kent & Falco 1988) and, combined with their 

observed position angles, the disk and bar can be weighed in- 
gravitational lensing is one of the only means by which these , 1 j_i /o 1 • ij_ 1 j_ T • ifir,o\ ivi- i 
° r, , , , r 1 , 7 7 dependently (bchmidt, Webster & Lewis 1998). Microlensmg 
properties ot high-redshitt galaxies can be probed. . , , 7 . , .,7 • .7 7 

7 7 r 7 n 01 the quasar images by compact objects withm the galaxy 

Conversely, lenses with low-redshift deflectors are also ^7 ^ r j_7 ■ /t ■ r 

permits the measurement ot their mean mass (Lewis & tr- 

win 1996; Wyithe, Webster & Turner 2000a), as well as a 
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ford 1992; Wyithe et al. 2000c), and even the transverse mo- 
tion of the galaxy relative to the line-of-sight to the quasar 
(Wyithe, Webster & Turner 1999). 

After Q 2237+0305, the second nearest deflector is the 
redshift 0.11 elliptical galaxy that lenses the radio source 
MG 1549+3047 (Lehar et al. 1993). Both the mass distri- 
bution and velocity dispersion of the galaxy have been mea- 
sured (Lehar et al. 1996), but the greater distance to the 
lens and less optimal source alignment place limits on the 
inferences that can be made from this system. 

It is clearly desirable to have a larger sample of lenses 
with low-redshift deflectors. Each would have the potential 
to provide as much insight as Q 2237+0305 and a sufficiently 
large sample could be used to extrapolate from the above re- 
sults to the galaxy and source populations (e.g. Keeton et 
al. 1998). The limitations of the various lens survey tech- 
niques are discussed in Section |^, with particular emphasis 
on the use of galaxy redshift survey (GRS) spectra as the 
primary data. Section ^ describes a simple model of the var- 
ious populations of objects that determine the viability of 
this method. In Section ^ specific predictions are made for 
various existing and planned redshift surveys, and the re- 
sults summarised in Section n. 



2 LENS SURVEYS 

Lensed quasars are very rare, with less than 50 known to 
date (e.g. Keeton & Kochanek 1996). At most 1 in ~ 100 
(and probably closer to 1 in ~1000) quasars brighter than 
m ~ 19 are lensed, which implies a surface density of 
^ 0.1 deg~^. This is several orders of magnitude lower 
than the integrated number counts of stars (~1000 deg"'^; 
e.g. Jones et al. 1991), galaxies (~ 100 deg~^; e.g. Mad- 
dox et al. 1990b), and unlensed quasars (~ 10 deg~^; e.g. 
Hewett, Foltz & Chaffee 1995) to the same limit. Thus lensed 
quasars represent between 0.001 per cent and 0.01 per cent of 
'bright' optical sources. Quasars (and hence lenses) represent 
a larger fraction of the radio source counts, and a number of 
large radio lens surveys exist (e.g. Patnaik et al. 1992; Myers 
et al. 1995). Unfortunately the lack of knowledge about the 
luminosity function and redshift distribution of the sources 
limits the power of any statistical results (Kochanek 1996b). 

Compounding these problems, considerable observa- 
tional effort is required to both identify candidate lenses 
and to then confirm their identity. Most are discovered in 
lens surveys (e.g. Burke, Lehar & Connor 1992; Crampton, 
McClure & Fletcher 1992; Maoz et al. 1992; Patnaik et al. 
1992; Maoz et al. 1993a,b; Surdej et al. 1993; Yee, Filip- 
penko & Tang 1993; Falco 1994; Kochanek, Falco & Schild 
1995; Jackson et al. 1995; Jaunsen et al. 1995) which are 
based on existing quasar (or extra-Galactic radio-source) 
samples. These catalogues in themselves represent a great 
deal of work, as initial multi-colour selection of candidates 
must then be followed by confirmation spectroscopy of each 
potential quasar. However there are now a number of large 
quasar samples in existence (e.g. Boyle, Shanks & Peterson 
1988; Warren, Hewett & Osmer 1991; Hewitt & Burbidge 
1993; Hewett et al. 1995) or in preparation (e.g. Loveday 
& Pier 1998; Boyle et al. 1999a,b), so any new lens survey 
can proceed from these object lists. Under this assumption, 
a lens survey then requires high-resolution imaging of be- 



tween 100 and 1000 fields (i.e. one for each quasar) per lens 
discovery. (Further observations of any probable lenses are 
usually required to establish the lensing interpretation be- 
yond doubt, but these represent only a small part of the 
overall task.) 

The efficiency of a lens survey can be increased by the 
use of morphological selection criteria, for instance choos- 
ing only highly-elliptical sources for further investigation. 
Kochanek (1991) found that, whilst the efficiency of such 
searches is quite high, they can be very incomplete, finding 
only ~ 20 per cent of all lenses. This conclusion was con- 
sistent with the results of the Automatic Plate Measuring 
(APM) lens survey (Webster, Hewett & Irwin 1988), which 
contains no lenses amongst ~ 2500 elliptical sources with 
quasar-like colours. 

Lens surveys also tend to miss some of the most inter- 
esting multiply-imaged quasars. These include those with 
high extinction and reddening from dust along the line-of- 
sight, and those with low-redshift lens galaxies. As discussed 
above, the latter can be very valuable indeed, but cannot be 
found in lens surveys they are not likely to enter the quasar 
surveys in the first place, due to both their morphologies (ap- 
pearing simply to be nearby galaxies at all but the highest 
resolutions) and colours (which are likely to be galaxy-like). 
Lenses with low-redshift deflectors can only be found by 
searching for quasar emission lines in the spectra obtained 
in GRSs. 

Indeed, Q 2237+0305 was discovered as a lens when it 
was noticed that one of the galaxies in the Center for Astro- 
physics (CfA) redshift survey (Geller & Huchra 1989) had 
an unusual spectrum. Follow up spectroscopy and imagin- 
ing confirmed that the object was in fact a lens. Whilst 
Q 2237+0305 represents, a posteriori, an unlikely align- 
menlj|], it is quite probable that other lenses will be discov- 
ered in redshift surveys, given that at least 1 in 10® redshifts 
will yield a lens (Kochanek 1992). Aside from the CfA sur- 
vey (which includes ~ 1.5 x 10^ objects), several other large 
GRSs have already been completed. The Las Campanas 
Redshift Survey (LCRS; Shectman et al. 1996) and the In- 
fraRed Astronomical Satellite (IRAS) 1.2 Jy survey (Fisher 
et al. 1995) have both measured upwards of 10'' galactic red- 
shifts, and so may contain a few lenses between them. With 
~ 9 X 10^ galaxies, the Muenster Redshift Project (MRSP; 
Schuecker, Ott & Seitter 1996) dwarfs the other completed 
surveys, but is probably less useful as a source of lenses, the 
redshifts having been measured from low-dispersion objec- 
tive prism spectra, as opposed to higher-resolution spectro- 
graphic data. In the future, the 2 degree Field (2dF) GRS 
(e.g. CoUess 1999; Folkes et al. 1999) wUl include ~ 2.5 x 10^ 
redshift measurements, and the Sloan Digital Sky Survey 
(SDSS; e.g. Szalay 1998; Loveday & Pier 1998) wiU con- 
tain about one million galaxies. The SDSS could yield 'sev- 
eral tens' of spectroscopically-discovered lenses (Kochanek 
1992), although it will also include high-resolution imaging, 
and thus should find a large number of lenses 'directly'. Im- 
portantly, these estimates are really only lower bounds on 
the number of lenses in GRSs, as previous calculations did 
not take into account two important observational effects. 

t Kochanek (1992) estimated the probability of a lens being dis- 
covered in the CfA survey was ~ 0.003, but see Section f^■l[ 
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The first is similar to the standard lensing magnification 
bias (e.g. Turner 1980), in which the number of lenses in a 
magnitude-limited sample is increased as intrinsically faint 
lensed sources are preferentially magnified into the survey. 
In a galaxy sample there is the additional effect that the 
magnified quasar images tend to be merged with the lens 
galaxy, making the composite object appear brighter. In this 
way (potential) lens galaxies up to one mag fainter than the 
survey limit can enter a GRS sample. As shown in Section ^ 
this effect can increase the number of lenses in a redshift 
survey by up to an order of magnitude. 

The second effect is complementary to this, as it only 
becomes important if the lens galaxy is much fainter than 
the quasar images (which is true of most lenses discovered 
to date). Such a lens would appear as several point sources, 
but might enter a GRS by virtue of appearing to be a single 
extended source in the low-resolution data from which sur- 
vey catalogues are generated. This is a potentially important 
possibility, but it is highly-dependent on the details of the 
survey selection effects and point-spread function (PSF) , as 
discussed in Sections ^ and ^. 

Assuming that a number of the objects in a GRS are 
lenses, and are hence observed spectroscopically as a mat- 
ter of course, a minimal number of confirmation observa- 
tions are required to find them. This method is only prac- 
tical if it is possible to differentiate between quasar spectra 
and galaxy spectra (as well as those of stars, white dwarfs, 
etc.), which, fortunately, is reasonably straight-forward, due 
to their prominent, broad emission lines. Other than the 
lenses, the only objects that will be selected in this way are 
any unlensed quasars that appear extended in the parent 
survey. These must also be re-imaged, and so the number of 
unlensed quasars in a GRS determines the efficiency of this 
type of lens search. 



3 SOURCE POPULATIONS 

The difficulty of finding gravitational lenses is principally 
a function of their rarity, so any assessment of search 
methods must include an estimate of the total number 
of sources which must be investigated per lens detection. 
The source popul atio ns considered are: stars (Section 3.1); 
galaxies (Section 3.2); unle nsed quasars (Section B.3); and 
lensed quasars (Section 3.4). Each population is defined by 



N{< m), the number of objects per unit solid angle brighter 
than magnitud^ m, which is shown for each of the four 
classes of objects in Fig. |l| 

In this simple model, the morphology of an object is 
considered to be completely defined by its ellipticity. 



+ My 



Myy)^ + AM^y 



M:CX + Myy ^J{Mxx - MyyY + 4Ma; 

where the moments, M, are given by 

Myy = {{ey - {6y}f} 



1/2 



, (1) 

(2) 
(3) 



■t The analysis presented here is reasonably generic, but in general 
Bj or B magnitudes are considered. 
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Figure 1. Cumulative number counts of stars (at the Galactic 
poles; dotted line), galaxies (dashed lines), unlensed quasars (up- 
per solid line) and lensed quasars (lower solid lines) as a function 
of Bj magnitude. The observed stellar number counts are taken 
from Stobie & Ishida (1987; filled squares) and Jones et al. (1991; 
empty squ ares ); the fit (from Bahcall & Soneira 1990) is described 
in Section [j.]] . The observed galaxy counts are from the APM 
galaxy survey (Maddox et al. 1990b; filled triangles) and Jones et 
al. ( 199 1; open triangles). The no-evolution galaxy models (Sec- 
tion [j.^ ) are given for Qmg = 1.0 and = 0.0; Quig = 0.3 and 
^Aq = 0-0; ^'^d f^niQ = 0.3 and = 0.7, although the number 
counts are almost independent of the cosmological model. The K- 
corrections of Kinney et al. (1996) are used for the lower triplet 
of dashed lines, whereas passive luminosity evolution is assumed 
to cancel the i^-corrections for the upper triplet of dashed lines. 
The observed number counts of unlensed quasars are taken from 
the compilation of data (the sum of the z < 2.2 and z > 2.2 sam- 
ples) presented by Hartwick & Schade (1990; filled circles), and 
the Fornax spectroscopic survey fP rinkwater et al. 1999; open 
circles); the fit given in Section |3.3t All lenses (i.e. irrespective 
of image separation or flux ratio) are included in the three lower 
solid lines, which assume: Qn 



1.0 and Qaq = 0-0 (lower line); 



= 0.3 and Qj^g = 
■ 0.7 (upper line). 



0.0 (middle line); and Qn 



0.3 and 



and 



M.y^({ey~{ey)){e,~(e,})}. (4) 

The angled brackets denote expectation values, defined by 



(h) 



h{e)f{e) Ae\ 



(5) 



where the integrals extend over the whole sky, and f{0) is 
the observed surface brightness of the object. The shape dis- 
tribution of each population is then defined by dj)/de|<m, 
the ellipticity distribution of all objects brighter than m. 
This is a significant simplification, the limitations of which 
are discussed where relevant. More complex methods (e.g. 
Jarvis & Tyson 1981; Godwin, Metcalfe & Peach 1983; Mad- 
dox et al. 1990a) provide more powerful methods of separat- 
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ing galaxies from stellar images, but are also more closely 
tied to individual surveys. The analysis presented here is 
reasonably general; Mortlock & Webster (2000c) present a 
more detailed and realistic simulation lensing in the 2dF 
GRS, which uses the sophisticated APM star-galaxy sepa- 
ration algorithm. 

3.1 Stars 

By far the most common bright point-like sources on the sky 
are stars, and their sheer numbers place limits on the effi- 
ciency of any search for lenses which includes point sources. 
Their density is not uniform across the sky, and thus it is 
preferable to perform any extra-Galactic survey away from 
the disk of the Milky Way, most commonly towards the 
Galactic poles. 



N 




0.2 0.4 0.6 0.8 1 



e 



3.1.1 Number counts 

The integrated number counts of stars in the B-band to- 
wards the Galactic poles can be approximated (to within 
~ 15 per cent) by 



iVs(< m) = (6) 

A^l lQ0.027{m-mgi ) 2^g0.083(m-ma2 ) 

[1 -I- 10-0.124{m-m,i)]3-l _|_ lQ-0.167(m-m,2 )] ^'^ ' 

where Ni ~ 3.9 x 10^ ~ 950 deg'^ m^i ~ 16.6, N2 ~ 
3.7 X 10'' ~ 910 deg~^ and ms2 18.0 (BahcaU & Soneira 
1980). This is shown as the dotted line in Fig. |l| along with 
the observed number counts (near the south Galactic pole) 
from Stobie & Ishida (1987) and Jones et al. (1991). 



3.1.2 Ellipticity distribution 

Despite the fact that all but the closest stars are un-resolved, 
most appear to have finite ellipticity due to noise effects. The 
APM survey is considered as a representative example of 
the low-resolution data used to generate galaxy survey cat- 
alogues, and so the stellar ellipticity distribution was chosen 
to match this, as described in Mortlock (1999). The func- 
tional form adopted was a two-dimensional Gaussian, given 
by 



dps 
de 



exp 



(e,)(<m) 



with the mean ellipticity over the range of interest (m ^ 22) 
well-approximated by 

1 /m-19^ 



(es)(< m) 



■ exp 



(^) 



(8) 



20 V 6 

This distribution is shown in Fig. ^ demonstrating both the 
variation with magnitude and the degree of overlap with the 
galactic ellipticity distribution (as given in Section 3.2). 



3.2 Galaxies 

The galaxy population not only dominates the number 
counts of extended sources, but also determines the number 
of lensed quasars, as the majority of deflectors are isolated 
galaxies. Three galaxy types are considered: spirals (S) and 



Figure 2. The ellipticity distributions of stars (and other point- 
sources) and galaxies, as labelled. The five distributions are for 
objects with m < 17 (soHd lines), m < 18 (dashed lines), m < 19 
(dot-dashed lines), m < 20 (dotted lines), and m < 21 (dot-dot- 
dot-dashed lines). The distribution of stellar ellipticities broadens 
and the galactic ellipticities becomes narrower towards fainter 
magnitudes, as the effects of noise begin to dominate. 



two classes of ellipticals (E and SO), with their relative num- 
bers given by vs — 0.69; ~ 0.12; and uso — 0.19 (Postman 
& Geller 1984). The co-moving galaxy luminosity function 
of each type is assumed to follow the same Schechter (1976) 
law, given by 

dng 
AM 



10 



■ ln(10)!^ng« 

-2/5 {M-M,){a + 1) 



exp [-10 



'2/5 (Af-Af,)] 



(9) 



where ng. = (0.016 ± 0.003)/i^ Mpc"^ gives the number 
density of all galaxies, M« = —19.4 ±0.1 in the B-band and 
a = -1.07 ± 0.05 (Efstathiou, Ellis & Peterson 1988). 

3.2.1 Number counts 

Assuming equation holds at all redshifts (i.e. no evolu- 
tion), the differential galaxy number counts for each popu- 
lation are given by 



diV, 
dm 



(7) 



AVo 
Az 



dn, 
dM 



Az, 



M=M(m,z) 



(10) 



where AVo/Az is the co- moving volume of a shell at redshift 
z (e.g. Carroll, Press & Turner 1992). Both the volumes and 
the distances are dependent on the current normalised den- 
sity of the universe, , and the normalised cosmological 
constant, Qaq. It is this dependence that allows both gravi- 
tational lensing statistics and number counts to be useful as 
cosmological probes, although the galactic number counts 
in the magnitude range considered here are almost indepen- 
dent of the cosmological model. The relationship between 
apparent and absolute magnitudes is 



M(m, z) = m — 25 — 5 log 



dL(0, z) 



IMpc 



K{z)-E{z), (11) 



where di^{0,z) is the luminosity distance (e.g. Carroll et al. 
1992), K{z) the /f-correction and E{z) the evolutionary cor- 
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rection, both of which are type-dependent. Two models are 
considered here: a non-evolving model with E{z) — and 
if-corrections as measured by Kinney et al. (1996); and one 
in which the AT-correction is cancelled out by the effects of 
passive luminosity evolution [i.e. E{z) — — Out to 
moderate redshifts the latter model gives a much better fit 
to the galaxy number counts (e.g. Fig. |l|; Maddox et al. 
1990b), and so is adopted as the default model. 

The cumulative number counts for each population are 
obtained by integrating equation (nol) to give 



iVg(< m) 



dVo dwg 
"d7 dM 



M — M(m',z) 

which is shown as the dashed lines in Fig. |l[ 



dzdm, (12) 



a 



d 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 














f s 

: --■'■■■.■/ ^ 

'-r^ 







0.1 



0.2 



0.3 



0.4 



0.5 



3.2.2 Ellipticity distribution 

Galaxies are the archetypal extended astronomical ob ject , 
but the ellipticity distribution of point sources (Section |3.l| ) 
is sufficiently broad that there is some overlap. This is espe- 
cially true for fainter objects - the ellipticity distributions of 
both stars and galaxies approach a noise-dominated distri- 
bution of intermediate mean ellipticity. A two-dimensional 
Gaussian is used to parameterise the distribution of galactic 
ellipticities [c.f. equation (M)], so that 



de 



exp 



(eg) 



where the mean ellipticity is given by 
- 20 



(eg) = 0.27 - 



25 



(13) 



(14) 



in the range of interest {15 ^ m ^ 23). This distribution 
is shown in Fig. |^ for a number of different magnitudes. 
For highly elliptical images the analytical function is too 
broad - e > 1 is permitted - but the behaviour of dpg/de 
is unimportant for e ^ 0.5, as galaxies are the only objects 
with such high ellipticities. 



3.2.3 Galaxy redshift surveys 

A generic GRS catalogue is assumed to consist of all images 
(in a certain region of the sky) which are brighter than the 
magnitude limit, miim, and have ellipticities greater than 
emin- For each source population the number over the whole 



l<"Hii 



which must then be 



sky is 7V(< mum)[l 
normalised by A'tot, the number of objects in the survey. 
The limitations of this model are discussed where relevant, 
but, despite its simplicity, it provides a reasonable estimate 
of the relative numbers of galaxies, stars and other objects 
that are observed spectroscopically in a galaxy survey. The 
ellipticity cut-off determines both the completeness, C, and 
efficiency, i5, of a GRS of a given magnitude limit; the resul- 
tant values of these parameters can then be used to assess 
the limitations of the model. 

A survey's completeness is the fraction of all galaxies 
which appear in the sample, and is simply 

C(emin,?niim) = 1 -Pg(< emin)|<mii„, (15) 

where pg(< e)|<m is given in equation (^|). Hence a survey is 
© 2000 RAS, MNRAS 000, 



Figure 3. The completeness, C, and the efficiency, E, of a GRS, 
as defined in equations and (^), shown as a function of the 
eUipticity cut-ofi^, e^i^. Both are shown for five different fimiting 
magnitudes: = 17 (solid lines); = 18 (dashed lines); 

'"lim = 19 (dot-dashed lines); nin^ = 20 (dotted fines); and 
miijn = 21 (dot-dot-dot-dashed fines). The galaxy number counts 
are calculated assuming that passive luminosity evolution cancels 
the /f-corrections, and an Einstein-de Sitter cosmological model 
is used. 



100 per cent complete if emin = 0, but becomes significantly 
incomplete for emin ^ 0.2, as shown in Fig. ^ 

The efficiency of a galaxy survey is defined as the frac- 
tion of objects in the survey which actually are galaxies, and 
so is dependent also on the stellar population. Using the 
stellar number counts and ellipticity distribution of equa- 
tions (Mj and (|7|), respectively, it is given by 



i?(emin,miim) 



(16) 



1 -Pg(< emin)|<mii„ 



Pg{< emin) t^n 



|<mii, 



iV3(<mii„) ' 



where Ng{< m) is given in equation (p^. This is shown 
for surveys of varying depth in Fig. ^ and the limitations 
of ellipticity as a means of star-galaxy separation are ex- 
posed. For a given magnitude limit the optimal choice of 
emin is close to the point where C = E, although slightly 
higher completeness may be preferable. Using emin — 0.15 
can result in completenesses and efficiencies of ~ 90 per cent 
for miim — 19, but more sophisticated techniques based on 
surface brightness profiles (e.g. Maddox et al. 1990a) can 
achieve ^ 95 per cent for both to fainter limits. 

The redshift coverage of a given GRS is of intrinsic in- 
terest, and it can also determine the lensing statistics within 
it (e.g. Kochanek 1992), depending on the selection effects. 
The redshift distribution of the galaxies is not independent 
of the ellipticity cut-off - the fainter and hence more distant 
galaxies are those likely to be omitted, reducing the high- 
redshift tail of the sample - but it is not a dominant effect, 
and is ignored here. From equations (^) and ([ic|), the a pri- 
ori probability that a galaxy of magnitude m has redshift z 
is 
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Az 



dVp dng 
dz dM 



M — M (m,z) 



r 

Jo 



dVp dng 
dz' dM 



Az' 



(17) 



M — A-/(m,z') 



A number of fitting formula liave been developed to approx- 
imate this function (e.g. Baugh & Efstathiou 1993; Brainerd, 
Blandford & Smail 1996), but the numerical integration is 
used here in order to look at the variation with cosmological 
model and evolution. The redshift distribution for the sur- 
vey is then obtained by integrating equation ( p^ over m and 
weighting by the differential number counts [equation (|l^)] , 
which gives 



Az 



r"i d}^ dng 
J_oo dz dM 



dm' 



M — M{m',z) 



roc rm dVp dn^ 
Jo J— oo dz' dM 



dm' dz' 



(18) 



Af = M(m',z') 



This is shown in Fig. ^ for both the E{z) = —K{z) model in 
(a) and the no-evolution model in (b). The mean redshift of 
the surveys increase with mum, but the influence of the K- 
corrections becomes increasingly dominant as the surveys 
go deeper. For ~ 20, the difference between the two 

models is approximately equivalent to a changing mum by 
one magnitude. 



3.3 Unlensed quasars 

The number counts of quasars brighter than the survey limit 
are important in this calculation as those with high elliptic- 
ities can contaminate any sample of candidate lenses. The 
faint end of the quasar luminosity function is important too, 
as it affects lensing statistics due to magnification bias. The 
redshift distribution of the sources is also necessary, as lens- 
ing probability increases so rapidly with source redshift (e.g. 
Turner, Ostrilcer & Gott 1984; Kochanek 1993a). 



3.3.1 Number counts 

The differential quasar number counts, as a function of mag- 
nitude, m, and redshift, z, are given by 



d^iVq _ 3[1 - iz/zm^f] 



qO 



AmAz 



6 + Zn 



■|^Qaq(m-mqo) _ J^Q.Sq (ni-mqo ) 



,(19) 



where mo — 19.0 ± 0.2 is the quasar break magnitude, Qq = 
0.9 ± 0.1 is the bright-end slope, /3q = 0.3 ± 0.1 is the faint- 
end slope, and iVqo = (4.1 ±0.4) x 10^ = 10±1 deg'^ (Boyle 
et al. 1988; Kochanek 1996a), and Zmax = 3.2 is chosen to 
match the redshift dependence of the Large Bright Quasar 
Survey (Hewett et al. 1995). Integrating equation (n9|) over 
redshift and magnitude gives the cumulative number counts 



Nq{< m) = 



r 

J — C 



qO 



■ Am, 



(20) 



which is plotted as the upper solid line in Fig. ^. 

The quasar and lens numbers calculated from these ex- 
pressions represent lower bounds on the true populations. 
Firstly, most quasar surveys are based on prior colour- 
selection (e.g. the UV-excess criteria of the Boyle et al. 1988 
sample), whereas any quasars found in redshift survey spec- 
tra are not subject to such cuts. There is also some uncer- 



tainty in the number of faint (m <; 21) quasars, as the deep- 
est quasar samples may be significantly incomplete. This 
possibility can be tested by 'complete' spectroscopic sur- 
veys, in which there is no morphological or chromatic pre- 
selection. Results from the deepest such sample, the Fornax 
spectroscopic survey (Drinkwater et al. 1999), are shown in 
Fig. |l|, but they do not differ greatly from the Hartwick & 
Schade (1990) compilation of data. Clearly a complete spec- 
troscopic survey to m ^ 21 would be of great interest. 

3.3.2 Ellipticity distribution 

As with stars, quasars can appear extended on photographic 
plates due to random noise. Their ellipticity distribution is 
likely to be similar to that of stars, but with more outliers, 
due to gravitational lensing and close quasar-galaxy associ- 
ations. As strongly-lensed quasars are considered explicitly 
here, the ellipticity distribution of unlensed quasars is as- 
sumed to be the same as that of stars, so that dpq/de|<m ~ 
dps/del<ni [equation (^ and Fig. |^. 



3.4 Lensed quasars 

The statistics of quasar lensing are most often framed in 
terms of the probability, pq, that a given quasar is lensed. 
Here, however, it is more useful to consider lensed quasars as 
a population of objects in their own right, distinct in partic- 
ular from the unlensed quasar population. The differential 
number of lenses with source redshift z and magnitude m is 
related to the lensing probability by 



dmdz 



Pi 



A'Nq 
AmAz ' 



(21) 



where A^ Nq/AmAz, is given in Section 3.3. The integral num- 
ber counts of lenses is then simply 



Ni{< m) = 



Am'Az 



Az dm'. 



(22) 



The first step in the calculation of pq is to find, Pq,g, the 
probability that the quasar (of magnitude m and redshift z) 
is lensed by a given galaxy. The lenses that are missed in 
follow-up observations (those with small image separations 
or high flux ratios) also have small ellipticities, and so are 
excluded from the model redshift surveys by the ellipticity 
cut described below and in Section ^. 

All three galaxy types are modelled as singular isother- 
mal spheres (Turner et al. 1984; Binney & Tremaine 1987; 
Kochanek 1994), which are consistent with lens statistics 
(e.g. Kochanek 1993a, 1996a) but cannot produce the more 
complex image configurations observed. Their mass distri- 
bution is characterised by a, the dark matter velocity dis- 
persion, which is marginally larger than the observed line-of- 
sight velocity dispersion, (T||, as shown by Kochanek (1994) 
and Mortlock & Webster (2000b). The line-of-sight disper- 
sion is related to M, the absolute magnitude of the lens 
galaxy, by the Faber- Jackson (1976) relationship for Es and 
SOs, and the TuUy-Fisher (1977) relationship^ for spirals. 
Both conversions have the same mathematical form, so 



3 The TuUy-Fisher (1977) relationship actually gives the circular 
speed, which is equal to \/2ct. 
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Figure 4. The redshift distribution of GRSs of varying magnitude limits: rrin^ = 17 (solid lines); miim = 18 (dashed lines); = 19 

(dot-dashed lines); = 20 (dotted lines); and = 21 (dot-dot-dot-dashed lines). In (a) the galaxies' /f-corrections are cancelled 

by passive luminosity evolution, whilst in (a) there is no luminosity evolution, and the i^T-corrections of Kinney et al. (1996) have been 
used. An Qmg = 1, ^Aq = cosmology is assumed in all cases as the cosmological model is unimportant at these depths. 



l.lcri 



-2/57 lM{m.g,Zg)-AU] 



(23) 



where 7e = 7so ~ 3.7 ± 1 and 73 = 2.6 ± 1 (de Vaucouleurs 
& Olson 1982) and a, = 225 ±20 km s"^ for E-type galaxies, 
= 205 ± 20 km s"^ for SO galaxies and a. = 143 ± 10 km 
s^^ for spirals (Efstathiou et al. 



Ignoring any selection effects for the moment, pq^g is 
given by the standard expression (e.g. Kochanek 1995; Mort- 
lock & Webster 2000b) 



where 



m'— m-t-5/2 log(/^tot) 



d^ito 



d2jVc, 



rfA(0,2) 



(24) 



(25) 



is the Einstein angle of the lens. Here dA(0, z) and dA(%, z) 
are the angular diameter distances from the observer to the 
quasar and the galaxy to the quasar, respectively (e.g. Car- 
roll et al. 1992; Schneider, Ehlers & Faico 1992). 

Integrating pq,g over the galaxy population gives 



Pq = 



dVb drig 
d^E dM 



Af — A/(mg,2g) 



Pq, 



r dirts: dZe 



(26) 



The choice of integration variables is arbitrary, but Zg and 
rrig are used here as they are observables. The resultant lens 
number counts are shown for three different cosmological 
models in Fig. |l|, but not all of these lenses are detectable 
in practice due to various observational limitations. 



3.4-1 Ellipticity distribution 

The vast majority of lensed quasars are effectively a small 
collection of stellar images, and whether they appear re- 
solved, extended, or point-like is strongly dependent on the 
PSF. The exception to this is when the lens galaxy is of com- 
parable brightness to the quasar images (e.g. Q 2237-f0305), 



in which case the source-deflector composite almost cer- 
tainly appears as an extended source; this case is dealt with 
separately. 

The singular isothermal sphere lens model, used above 
to calculate the number of lenses, is too simple to reproduce 
the observed variety of image configurations. A quadrupole 
term - either intrinsic ellipticity of the lens or external shear, 
as used here - is required to generate the four- (and five-) 
image lenses observed. The possibility of a small core ra- 
dius is less important, as in most cases it only results in an 
(additional) demagnified central image. In the presence of a 

dimensionless external sheaip] 70, but with no core radius, 
the lens equation is (e.g. Kovner 1987; Schneider et al. 1992) 



1-70-^ 



(1 



(27) 



where /3 = {(3x,(3y) is the angular source position, 6 = 
{dx, Oy) is the angular image position. This can be solved us- 
ing the method outlined by Schneider et al. (1992), in which 
the radial components of the image positions are given by 
the roots of a one-dimensional equation. 

In order to obtain the ellipticity distribution of lensed 
quasars, Monte Carlo simulations of ~ 10^ lenses were gen- 
erated for several different deflector models. The integra- 
tion over the quasar population was performed by using pq 
[equation (^] and d'^ Nq/dzdm to draw sources randomly 
from population of lensed quasars, and a lens galaxy was 
similarly-generated for each source. The source position, (3, 
was biased in favour of close alignments using the magni- 
fication bias described above, but no other selection effects 
were included. For each source-lens pair the lens equation 
was solved to give image positions, Oi and magnifications. 



^ The use of an elliptical lens model is more realistic, but more 
expensive computationally. A shear of magnitude 70 is roughly 
equivalent to an ellipticity of 870 (e.g. Keeton & Kochanek 1996). 
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Figure 5. Dimensionless forms of the PSF, /scc(t*), under the 
assumption of a Moffat (1969) profile. The profile is shown for 
different values of the rj-parameter used in the definition of the 
PSF [equation (^)]: r] = 2.0 (solid line); rj = 2.5 (dashed line); 
rj = 3.0 (dot-dashed line); rj = 5.0 (dotted line); and ?} — > oo (i.e. 
Gaussian seeing; dot-dot-dot-dashed line). 



Hi, with 1 < i < Ni, where Ni is the number of images 
(either two or four in this model). 

The appearance of a given lens in the survey data is 
determined by the PSF, denoted by /sec (6) ~ fscc{9), which 
is normalised so that JJ fsccid) dO^ = 1. The contribution 
from the galaxy is usually negligible (See below.), and sur- 
face brightness-related selection effects are also ignored here 
(but treated explicitly in Mortlock & Webster 2000c). The 
surface brightness of a lens is then 



(|6I-6I,;| 



(28) 



where F is the total observed flux of the quasar images. From 
this the moments and ellipticity [as defined in equation (|l|)] 
can be found. 

In a similar Monte Carlo simulation, Kochanek (1991) 
assumed Gaussian seeing, ignoring the extended wings of the 
PSF (e.g. King 1971). A more general Moffat (1969) profile 
is used here, for which the PSF is 



(r,~l)(2i/"-l) 



which can be integrated to give 



(2'/" - 1) 



Fscci< 9)^1- 



l + (2 



1) 



.(29) 



(30) 



where 6^5 is the full width at half-maximum of the seeing disc 
and rj determines its shape, as shown in Fig. ^. The profile 
cannot be normalised if »7 < 1, and in the limiting case of 
77 — > oo it becomes a Gaussian. The default value assumed 
by the IMAGE reduction and analysis facility software 
(Tody 1986) is rj = 2.5, but results such as those of Saglia 
et al. (1993) imply that ?7 ~ 3. Conversely the areal profiles 
of stellar images in the APM survey shown in Maddox et al. 
(1990a) are close to log- linear, suggesting rj is very high. 

Using this PSF the ellipticity can be found analytically, 
although the result is rather cumbersome, except if Ni = 2. 



In this case it can be assumed, without loss of generality, 
that both images lie on the x-axis, and hence equation (uh 
reduces to 



1 



{{ey - {ey)f) 



-I 1/2 



Using equation (|2s|), this becomes 

e = 



1 - 



l + 2(r/-2)(2i/''-l) 



10 



2Am/5 



(1 + 102^^ 



m/5\2 



A9_ 

97/2 



(31) 



(32) 
1 1/2 



where is the angular separation of the images and 
Am is their magnitude difference. This is only defined for 
rj > 2, and reduces to the simpler Gaussian form found by 
Kochanek (1991) in the limit of 77 — » 00. Fig. ^ shows e as 
a function of A6 for several values of Am and two extreme 
values of ry. For a given image separation and seeing, e de- 
creases with increasing Am, the pair becoming more like a 
single point source. The ellipticity also decreases with de- 
creasing ?7, the wings of the PSF 'diluting' the ellipticity of 
the two images. More generally, it is clear that small sep- 
aration {A6 ^ 1 arcsec) lenses stand little chance of being 
classified as anything but stars. 

Using the above formulation to calculate e for each sim- 
ulated lens gives a distribution, dpi/de|<m, that is depen- 
dent mainly on 70 and the PSF. There is only weak depen- 
dence on the cosmological model, the magnitude limit of the 
survey, and the scale velocity dispersion of the galaxies (as 
(Ae) (X al; Kochanek 1993b). Fig. g shows the cumulative 
ellipticity distribution, pi(< e)|<m, for various values of the 
model parameters. The critical value is 1 — pi(< eniin)|<m, 
as this is the fraction of lenses that would enter any given 
GRS. From the completeness and efficiency arguments in 
Section 3.2, 0.1 ^ fimin ^ 0.2; this range is delineated by the 
vertical dashed lines. An external shear of 0.2 can reduce the 
number of 'elliptical' (i.e. e > emin) lenses by up to a factor 
of four, but the effects of the PSF are clearly more impor- 
tant. Increasing 9s from 2 arcsec to 3 arcsec decreases the 
number of elliptical lenses by ~ 30 per cent in Gaussian see- 
ing, but results in an order of magnitude decrease if ~ 2.5. 
With 77 ^ 3 the PSF is broader, and 0.1 <, pi{< e)|<™ <, 0.3 
(for 2 arcsec seeing) or 0.05 ^ Pi{< e)|<m ^ 0.2 (for 3 arc- 
sec seeing) . The strong dependence on the shape of the PSF 
shows that the wings do contribute significantly to the el- 
lipticity of multiple images, contrary to the assumption of 
Kochanek (1991). The variation in PSF between different 
sets of observational parameters thus places limits on the 
accuracy of any generic predictions of the number of lenses 
that will enter a GRS. 



3.4- S Lenses with visible deflectors 

In the most pessimistic scenario presented above, only about 
5 per cent of all lensed quasars are sufficiently elliptical to 
enter a GRS as candidate galaxies, in which case there would 
be only limited use in searching GRS spectra for potential 
lenses. However the simulation did not take into account the 
light from the deflector. If the lens galaxy is bright enough 
for the source-deflector composite image to appear extended. 
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Figure 6. The ellipticity, e, of an image pair witii angular separation A8 (with 0s being the seeing). Results are shown for various values 
of the magnitude difference between the images: Am = (solid lines); Am = 1 (dashed lines); Am = 2 (dot-dashed lines); Am = 3 
(dotted lines); and Am = 4 (dot-dot-dot-dashed lines). A Moffat (1969) PSF is assumed, with r] = 2.5 in (a) and jj — > oo (i.e. Gaussian 
seeing) in (b). The horizontal dashed lines are at e = 0.1 and e = 0.2, reasonable limits for the ellipticity cut-off used to select objects 
for galaxy surveys. 




Figure 7. The cumulative distribution of ellipticities of lensed quasars, pi(< e)|<m, in seeing of 6s = 2 arcsec (a) and 9b = 3 arcsec (b). 
In each case an Einstein-de Sitter cosmology is assumed and the magnitude limit is m = 19. The lens model is a singular isothermal 
sphere with an external shear of: 70 = 0.0 (solid lines); 70 = 0.1 (dashed lines); and 70 = 0.2 (dot-dashed lines). The three lines for 
each lens model are for different values of the 7;-parameter used in the definition of the PSF [equation (^9[)]: rj = 2.5 (upper lines); rj = 3 
(middle lines); and »; — > 00 (i.e. Gaussian seeing; lower lines). The vertical dashed lines indicate reasonable values for the ellipticity cut-off 
used to select objects for a galaxy survey, and this can be used to obtain the fraction of lenses retained. 



but the quasar is also bright enough that the composite spec- 
trum is distinguishable from that of an isolated galaxy, the 
lens will both enter the redshift survey and be detectable 
spectroscopically. Kochanek (1992) performed an 'initial, or- 
der of magnitude' estimate of the statistics of such lenses, 
but the contribution of the quasar light to the observed mag- 
nitude of the compound object was neglected. This is an 
important effect, as it increases the effective depth of the 
lens survey by ~ 1 mag, and can result in up to an order of 
magnitude more lenses. 

Defining mg as the magnitude of (just) the galaxy, and 
JTiq as the total, lensed magnitude of all the quasar images, 
the observed magnitude of a lens is given by 



m = -5/2 log ^ ^Q-2/5m^^ ^ (33^ 

As illustrated in Fig. ^, a lens with these properties will be 
included in a redshift survey of limiting magnitude if 
the following three criteria are met: 

(i) The lens must be bright enough to enter the GRS at 
all, i.e. m < mum, which is shown by the curved line in 
Fig. ^ This gives the extra depth relative to the redshift 
survey proper, for which the Hux cut-off is rUg < mu^. 

(ii) The quasar must be sufficiently bright, relative to the 
lens galaxy, that its emission lines are detectable in the GRS 
spectrum. Following Kochanek (1992), this requirement is 
modelled as mg — rriq < Am-qg, which is given by the upper 
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of the two diagonal cuts in Fig. ^ Due to the strong, broad 
emission lines present in most quasar spectra (e.g. Peterson 
1997), Arriqg ~ 2 is assumed unless stated otherwise. This 
value is also supported by the discovery of Q 2237+0305 as 
described in Kochanek (1992). 

(iii) The galaxy must be sufficiently bright that the the 
lens is classified as a galaxy in the low-resolution data from 
which the GRS is selected. This is assumed to be the case 
if JTiq — rrig < Arrigq, which is shown as the lower diago- 
nal cut in Fig. ^. The value of Amgq is somewhat uncer- 
tain, and could even be negative^ if galaxies with super- 
imposed stellar components are reliably removed from the 
sample. In the case of Q 2237-1-0305, rriq ~ mg + 1, but the 
galaxy appears unremarkable at low resolution, suggesting 
that Arrigq ~ 0. The simulations of the APM galaxy survey 
presented in Mortlock & Webster (2000c) further imply that 
Arrigq 1. In other words, the galaxy need contribute 
only half the flux of a lens for it to be classified as non-stellar, 
and hence included in a GRS. If the 'elliptical' lenses consid- 
ered previously enter GRSs, they could be included in this 
formalism by taking Arrigq oo. 

It is now possible to calculate Pq,g, the probability that 
a given galaxy (with redshift Zg and magnitude rrig) lenses 
a given quasar [with redshift z, and lensed magnitude rriq, 
defined in terms of m and rrig, in equation (^^]. Adjusting 
equation (pll), the lensing probability becomes 



8 d^iVq 



4 L M?ot d^dm' 



m' = --5/2 log(lO-2/5 r 



-2/5 1 



(34) 
d/itot , 



if rrig^min < rrig < rrig. max, but is zero if nig is outside these 



limits. These, in turn, are given by 



ain = ?ri + -log(l + 10 



-2/5 Amqg 



(35) 



and 



= m+|log(l + 10'/^^'"-) 



(36) 



The integral is along a line of constant m, as shown in Fig. 
Integrating pq,g over the galaxy population yields 




Figure 8. The range of deflector and source magnitudes (rrig 
and rriq, respectively) probed by a GRS of magnitude limit mw^. 
Lenses with —5/2 log(10~-^''^ ™<i -I- lO"^''^"'^) > are omit- 

ted as the combined flux from both the quasar and lens galaxy is 
below the flux limit. The survey is assumed to have Arriqg = 2 (all 
lenses with rriq — rrig > Arriqg being omitted as the quasar is too 
faint to be detectable spectroscopically) and Arrigq = 1 (all lenses 
with rrig — riiq > Arrigq being omitted as the galaxy is too faint to 
make the lens appear extended) . Hence only those lenses which lie 
in the un-hatched region are included. The heavy solid lines show 
the boundaries of the equivalent region used by Kochanek (1992) - 
the inclusion of the lenses for which nig > miim is far more impor- 
tant than exclusion of the lenses in which the galaxy is too faint. 
This last point is illustrated by the scatter plot within the al- 
lowed region, which shows the distribution of galaxy-quasar pairs 
generated from their respective number counts [equations ( [l2[ ) 
and (p^)] . For the scatter plot mn^^ = 19 was assumed, and the 
quasars, being strongly lensed, are assumed to be magnified by a 
factor of 4 on average. 



d^iVi 
dmdz 



(37) 



dVb drig 
dZ dM 



M=M{mg,Zg) 



Pq.gdrrig dzg, 



where the integral over the defiector magnitude is restricted 
again by the requirement that the overall magnitude of the 
lens be m. Finally the cumulative number counts of such 
lenses are 



Ni{< m) 



dmdz 



dz. 



The values of A^'i implied are discussed in Section |4.1 
is immediatel 
pairs in Fig. ! 



(38) 



but it 



clear from the distribution of galaxy-quasar 
that the majority of the lenses are included 



-Arriqg, then the only lenses included would be 
those with quasars too faint to be detectable, and the lensing 
probability would be zero. 



only if both the quasar and galaxy fluxes are explicitly ac- 
counted for in the calculation. 



4 RESULTS 

Combining the various populations of objects considered 
in Section it is possible to make predictions about the 
number of lenses found in particular redshift surveys (Sec- 
tion [l.l| ) as well as the usefulness of this search method 



(Sections h.i 4.3 and 4.4) 



4.1 Number of lenses 



From the number counts given in Section |^, the number of 
objects of type t which enter a given redshift survey is 



(39) 



Nt{< miim)[l -Pt(< Cmin) 

Cmin) I <?Ti^jj.j^] 
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where A'tot is the total number of objects in the survey, 
miim the magnitude hmit and emin the eUipticity cut-off. 
The sum is over all the populations (i.e. t — s, g, qor 1, for 
stars, galaxies, unlensed quasars and lenses, respecitvely), 
although only stars and galaxies make a significant contri- 
bution to the denominator. 

Fig. 1^ shows the number of 'elliptical' lenses (and other 
objects) per 10* redshifts as a function of (chosen to 

facilitate comparison with Kochanek 1992) for a number of 
different observing conditions. The lens fraction remains es- 
sentially constant at about 1 in 10^ (a) or 1 in 10* (b), de- 
creasing very slowly as mum increases, due to the decreased 
magnification bias. For bright surveys this is far greater than 
the numbers of lenses with bright defiectors considered by 
Kochanek (1992), but it is not at all certain that the 'ellip- 
tical' lenses will enter a typical GRS. 

Multiply-imaged quasars with lens galaxies close to the 
survey limit should definitely enter galaxy samples, and so 
they give a lower bound on A^i. As can be seen from Fig. |lo[ 
these lenses are considerably rarer than the 'elliptical' lenses 
at bright magnitudes, but their numbers increase rapidly 
with miiin, due t o t he steepness of the quasar luminosity 
function (Section 3.3). Fig. hu can be compared directly with 
the results of Kochanek (1992), although Ni is higher by a 
factor of ~ 5, due to the inclusion of the quasars' light in 
the calculation. This immediately implies that the majority 
of any such lensed quasars found in redshift surveys should 
have lens galaxies which are actually fainter than the survey 
limit, as can be seen from Fig. |^. Other than the depth of the 
survey, M is determined mainly by Arriqg, with Amgq some- 
what less important. The strong dependence of A^i on Amqg 
is again purely a function of the quasar number counts, but, 
being a measure of how sensitive a redshift survey is to the 
presence of quasar emission features, it can be determined 
'experimentally' by running the survey software on the spec- 
tra of simulated lenses. The value of Amgq is more difficult 
to determine; fortunately it does not greatly affect the lens 



Of the other large redshift surveys that have already 
been completed, both the European Southern Observatory 
Slice Project (ESP; Vettolani et al. 1997) and the Cana- 
dian Network for Observational Cosmology (CNOC2) sam- 
ple (Yee, EUingson & Carlberg 1996) are too small to present 
any real likelihood of lensing. The million galaxies in the 
MRSP (Schuecker et al. 1996) could include a number of 
lenses, although only objective prism spectra were used 
for the redshift determination. The selection of candidate 
quasars from low-grade spectra is an established technique 
(e.g. Hewett et al. 1995), but this would represent a consider- 
able undertaking, as the low resolution data would generate 
a large number of false candidates lenses for which further 
observations would be required. Thus the most promising of 
the completed galaxy surveys is the LCRS (Shectman et al. 
1996) - it is larger than the CfA survey, as well as being 
deeper ((m) ~ 19), and should contain one or two lenses 
(c.f. Fig. @. 

The two surveys listed in Table ^ which have not been 
completed are also the most ambitious, and are likely to yield 
the most lenses. The 2dF GRS (e.g. CoUess 1999; Folkes et 
al. 1999) should contain at least 10 lenses, and as many as 50 
if the cosmological model and observational parameters are 
optimal. Either way, it would represent the largest sample of 
lensed quasars generated from a single survey. However the 
2dF instrument has very small optical fibres (effective radius 
1 arcsec), and is unusually sensitive to a number of surface 
brightness-related selection effects. The impact of these on 
Ni is studied in more detail in Mortlock & Webster (2000c). 

The SDSS (e.g. Szalay 1998; Loveday & Pier 1998) 
is four times the size of 2dF, and the spectra will be of 
considerably better quality. Assuming Arriqg ~ 4 implies 
that Ni J; 100. Further, the survey will also include high- 
resolution imaging of ~ 10* deg^, which should allow the 
morphological identification of an even larger number of lens 
candidates. The number counts shown in Fig. ^ then imply 
that the SDSS should contain between 200 (if flmp — 1.0 



statistics, provided that Arriqg ^ 1. The cosmological model and Qaq = 0.0) and 1000 (if Q^o = 0.7 and = 0.3). 
also has little impact on the number of lenses with visi- 
ble deflectors, as the lens galaxies are so close-by that the , ^ , 

, , „ , ,. , 4.2 Control samples 

observer-source and deflector-source distances are approxi- 

mately equal (Kochanek 1992). 

These above results are summarised for a number of 
real redshift surveys in Table 0. None of the surveys have 
been systematically searched for lenses, although all the CfA 
survey spectra were examined by eye, and it was process 
that led to the discovery of Q 2237-1-0305 as a lens. If 'el- 
liptical' lenses were included in the survey, the detection 
of a lens is not unexpected; however the CfA GRS is so 
bright {{m) ~ 15) that multiple point sources should have 
been reasonably easy to remove from the sample. Assuming 
only those lenses with visible, low-redshift deflectors were 
included, the expected number of lenses in the CfA survey 
is ~ 0.03. In other words, about 1 in 30 CfA-like galaxy 
samples should contain a spectroscopic quasar lens. This is 
considerably higher than the previous estimates of ~ 0.0002 
(Huchraet al. 1985) and ~ 0.003 (Kochanek et al. 1992); the 
reason for this increase is simply the inclusion of the quasars' 
light in the calculation. Note also that Q 2237-1-0305 is ac- 
tually fainter than the survey limit, which is consistent with 
the conclusion that galaxies with rrig ^ rriiim dominate the 
lens statistics. 



Kochanek (1992) suggested two control samples - fore- 
ground stars and unlensed quasars - from which the effi- 
ciency and completeness of the lens search could be cal- 
ibrated. However these objects are only useful controls if 
they enter the GRS mainly through chance superpositions 
with survey galaxies. Unfortunately, most of the stars and 
unlensed quasars in the survey are not there due to chance 
alignments, but misclassification (e.g. Fig. For instance, 
the limitations APM star-galaxy separation algorithm (e.g. 
Maddox et al. 1990a) imply that, of the ~ 2.5 x 10^ objects 
in the 2dF GRS, several thousand will be stars, compared 
to only several tens of star-galaxy superpositions. With the 
exception of the brightest surveys (niiim ^ 18), the same is 
also true for unlensed quasars. 

This effectively leaves this lens search technique without 
any simple calibration, which is not important if the aim is 
simply to discover new lenses, but must be addressed if lens- 
ing probabilities are to be determined. One possibility is to 
perform both morphological and spectroscopic calibration a 
posteriori by simulating quasar-galaxy lenses and analysing 
them using the software actually used for the survey. This 
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Figure 9. The expected number, A^t, of galaxies, stars, unlensed quasars and lensed quasars in a GRS as a function of the magnitude 
limit, miinj. The number of lenses is calculated using the spherical lens model (i.e. 70 = 0) in three cosmologies: = 1-0 and = 0.0 
(lower solid line); = 0.3 and = 0.0 (middle solid line); and Omg = 0.3 and = 0.7 (upper solid line). The GRS is selected 

using an ellipticity cut emin = 0.15, and normalised to A'^tot = 10**. The parent survey is assumed to have 3 arcsec seeing and a Moffat 
(1969) PSF with r\ = 2.5 in (a) and — > 00 (i.e. Gaussian seeing) in (b). 
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Figure 10. The expected number of lensed quasars 
magnitude limit, 
and -5/21og(10-2/5 



A^i, per 10^ measured redshifts, shown (as in Kochanek 1992) as a function of 



+ 10-2/5 mg) ^ 



are included. (See Section 3.4 for a full explanation of these selection criteria.) In each 

(lower lines); Amqg = 2 (middle 



panel the three sets of lines represent different values of the spectral prominence of quasars: Arrtqg ■ 
lines); and Amqg = 4 (upper lines). The morphological threshold for galaxy selection is assumed to be Arrigq = in (a) and Arrigq = 1 
in (b). For each set of observational limits, three cosmological models are shown: Omo = 1-0 and = 0.0 (solid lines); Omo = 0.3 and 
= 0.0 (dashed lines); and = 0.3 and f2yvg = 0.7 (dot-dashed lines). 



also has the advantage of taking into account all of the po- 
tential biases in the survey without the need for explicit 
modelling. 



4.3 Efflciency and completeness 

The sheer number of quasars that must be re-imaged per 
lens discovery in a conventional lens survey is prohibitive 
(e.g. Section y), and it is one of the major reasons that so 
few lenses are known. The use of morphological selection, 
implicit in a GRS, can increase the efficiency by removing 
most of the unlensed quasars from the survey, but there is 



also a significant reduction in the completeness of the lens 
sample. 

The completeness is simply C — 1 — pi(< emiii)l<mii„ij 
where emin is the ellipticity cut of the survey and its 
magnitude limit. As shown in Fig. ^ this is quite variable, 
with C ranging between 0.05 and 0.5, depending mainly on 
the form of the PSF. However a low value of C is not neces- 
sarily a problem, provided that its value is reasonably well 
known. Note that C is unambiguously defined only for the 
'elliptical' lens sample; the lenses with rrig ~ mq are not 
drawn from a well-defined parent population. 

More important is the efficiency, which is defined as 
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Table 1. The frequency of lensing in existing and planned GRSs. 



Survey 


Reference 










A^tot 




Ni 






M.min 


















= 1 




= 0.3 


^mg — 1 


















= 




= 0.7 


f^Ao = 
















Tj = 2.5 




= 2.5 






CN0C2 


Yee et al. (1996) 


R = 


21.5 


Bj 


= 22.8 


2.6 X 10^ 


0.01 


0.03 


0.08 


0.2 


0.3 


ESP 


Vettolani et al. (1997) 


Bj = 


19.4 


Bj 


= 19.4 


3.3 X 10^ 


0.04 


0.1 


0.2 


0.7 


0.1 


CfA 


Geller & Huchra (1989) 


-Bzwicky = 


15.5 


Bj 


= 15.5 


1.5 X 10* 


0.1 


0.4 


0.8 


2.3 


0.03 


LCRS 


Shectman et al. (1996) 


R = 


17.7 


Bj 


= 19.2 


2.6 X 10* 


0.3 


0.8 


2.3 


5.7 


0.9 


2dF 


CoUess (1999) 


Bj = 


19.45 


Bj 


= 19.45 


2.5 X 10^ 


2.5 


7.5 


18 


50 


10 


MRSP 


Schuecker et al. (1996) 


J 


= 20 


Bj 


= 19.5 


9.0 X 10^ 


4.5 


18 


36 


90 


27 


SDSS 


Szalay (1998) 


'''Gunn = 


18.2 


Bj 


= 18.9 


1.0 X lO'* 


10 


40 


90 


250 


100 



A summary of the number of lenses expected in several completed and future redshift surveys. For each survey: is the magnitude 

limit in the band indicated (with approximate conversions to Bj, based on the redshift coverage of the survey); A'^tot is the total 
number of redshifts measured; A^i is the expected number of lenses; and A^i is the expected number of lenses if only those with 
visible detectors are included in the survey. The variation of Ni with cosmological model and PSF is given explicitly, but the spherical 
lens model and 3 arcsec seeing are assumed. On the other hand, A^i^min is strongly dependent on neither the cosmology, nor the PSF 
in the relevant magnitude range, and so is only shown for the Einstein-de Sitter model. In all cases Ni is calculated assuming 
Amgq = and Arrtqg = 2, except for the MRSP (which uses objective prism spectra, and so has Amqg ~ 0) and the SDSS (which 
will obtain high-quality spectra with Amqg ~ 4). 



the number of high-resolution re-observations required per 
lens discovery. For a redshift survey-based lens search it is 
[c.f. equation (|l^)] 

A''l(< miim)pi(< emin) j<m,i„ 



-Egrs = 



A'"q(< miim)pq(< emin) \<r. 



(40) 



as only those objects with quasar-like spectra need to be re- 
imaged. As shown in Fig. h]], this is generally much higher 
than the efficiency of a conventional lens survey. 



Nii< mii„ 



lens survey — Ar / \ ' \ J 

iVq(< miim) 

where it is optimisticly assumed that the lens sample is 
essentially complete. For bright surveys -Eqrs ~ 1 as the 
star-gala:xy separation techniques are so reliable, although 
close quasar-galaxy associations have not been included. For 
fainter magnitude limits E decreases for both search meth- 
ods, due to the lower lensing fraction, and for mum 20 the 
efficiency of redshift surveys approaches that of lens surveys. 

If only lenses with bright deflectors enter a galaxy sur- 
vey, then iJoRS is reduced by about an order of magnitude, 
but is still greater than Sicns survey for miim ^ 20. Hence the 
efficiency of lens searches based on GRS spectra is almost 
always greater than that of conventional lens surveys. 

4.4 Deflector redshift distribution 



By evaluating the lensing probability (Section i.4) for a fixed 
deflector redshift, the distribution of lens galaxy redshifts, 
dp/dzg, can be found; this is shown in Fig. [13. The 'elliptical' 
lenses have the broad range of deflector redshifts expected of 
a conventional lens survey (e.g. Turner et al. 1984), and, in 
particular, the expected number of nearby deflectors is very 
small. The distribution is very different if only lenses with 
visible lens galaxies are selected - it roughly matches the 
overall redshift distribution of a GRS with a limiting mag- 
nitude of miim + Arriqg (c.f. Fig. Hence the 2dF redshift 
survey should yield one or two lenses with % ^ 0.1, and 
the SDSS could contain about ten lenses with such nearby 
deflectors. 




Figure 11. The detection efficiency of a conventional lens survey, 
^lens survey (lower lines), and that of a lens search based on a 
redshift survey, -Eqrs (upper lines), as a function of the limi ting 
magnitude of the survey, mn^. Both are defined in Section [4.3| , 
and -EgrS is shown for two values of r] used in the definition of 
the Moffat (1969) PSF: r] = 2.5 and 77 — > oo (Gaussian seeing), as 
labelled. In each case three cosmological models are used: f2mo = 
1.0 and = 0.0 (solid lines); = 0.3 and = 0.0 (dashed 
lines); and fluiQ = 0.3 and ^An = 0.7 (dot-dashed lines). 



5 CONCLUSIONS 

Due to the imperfect discrimination between galaxies and 
other celestial objects, gravitationally-lensed quasars enter 
GRS catalogues, and should then be detectable spectroscop- 
ically. The gravitational lens Q 2237-1-0305 was discovered 
in this way (Huchra et al. 1985), and Kochanek (1992) pre- 
dicted approximately one lens per 10*' redshifts measured. 
However the inclusion of the quasars' light in the calcula- 
tion of the lens galaxies' magnitudes increases the number 
of lenses by up to an order of magnitude, as illustrated in 
Fig. 1^. Another possibility is that many lensed quasars enter 
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Figure 12. The normalised distribution of deflector redshifts, Zg, 
of lensed quasars discovered in GRSs. For the 'elliptical' lenses 
(broader curves) the distribution is determined solely by the cos- 
mological model, whereas the extent of the distribution for the 
lenses with bright deflectors increases with the depth of the sur- 
vey. The two sets of curves for each type of lens show the dis- 
tributions with miinj = 18 and = 20. For each survey 
Amqg = 2, and results are presented for three cosmological mod- 
els: HniQ = 1.0 and Oaq = 0.0 (solid lines); Hmo = 0.3 and 
f^Ao = 0.0 (dashed lines); and Cmo = 0.3 and SIaq = 0.7 (dot- 
dashed Hues). 

GRSs as the quasar images combine to have a high elliptic- 
ity, but are also unresolved in the low-resolution (e.g. plate) 
data from which candidate galaxies are selected. If a signifi- 
cant fraction of such lenses are observed spectroscopically in 
GRSs, more new lenses will be discovered in redshift surveys 
than are known to date. 

The number of lenses expected in various existing and 
planned surveys is given in Table 0. It can be seen that A'^i 
varies greatly with both the form of the PSF of the par- 
ent survey and the cosmological model (in the case of the 
'elliptical' lenses), which places fundamental limits on the 
accuracy of these predictions. For instance, the current gen- 
eration of surveys (with up to ~ 2 x 10^ galaxies) should 
contain several lenses between them, but the above uncer- 
tainties, combined with shot noise, make this a rather weak 
prediction. Looking ahead, the 2dF GRS (with ~ 2.5 x 10^ 
galaxies to a limit of ms, = 19.5) is already well under- 
way, and should contain at least 10 new lenses, and up to 
50 if observational conditions are favourable. Howver, it is 
particularly sensitive to surface brightness-related selection 
effects, and a more detailed simulation of lensing in the 2dF 
redshift survey is presented in Mortlock & Webster (2000c) . 
Finally, the SDSS can be expected to contain over 100 spec- 
troscopic lensed quasars, along with even more discovered 
by more conventional methods. 

Many, and possibly all, of the lenses discovered spectro- 
scopically in redshift surveys, will have deflector galaxies at 
low redshifts, comparable to the depth of the survey proper. 
These have the potential to be the most important products 
of such a lens search, as the proximity of the deflector, com- 
bined with the information provided by the lensing event, 
can provide a number of unique insights, as exemplified by 
Q 2237+0305. The 2dF GRS should contain several low- 



redshift lenses, and the SDSS about 10. Of course not all 
of these will have Q 2237-|-0305's wonderful combination of 
source and defiector properties, but the possibility of even 
one similar system is tantalising, to say the least. 
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